Density fluctuations of polymers in disordered media 
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We study self avoiding random walks in an environment where sites are excluded randomly, in 
two and three dimensions. For a single polymer chain, we study the statistics of the time averaged 
monomer density and show that these are well described by multifractal statistics. This is true even 
far from the percolation transition of the disordered medium. We investigate solutions of chains in 
a disordered environment and show that the statistics cease to be multifractal beyond the screening 
length of the solution. 



INTRODUCTION 



In many experimental situations, linear polymers are 
present in a disordered medium. For example, a free 
polymer chain inside an elastomeric network, or a DNA 
molecule in an agarose gel. The properties of single 
chains in such situations have been the subject of numer- 
ous theoretical and numerical work From a theoret- 
ical perspective, the problem is of great interest in part 
because of its strong connection to the phenomenon of 
quantum mechanical localization [Toj . The greens func- 
tion for a single ideal polymer chain in a potential is 
mathematically identical to that of an electron in the 
same potential, but in imaginary time (Tlj . Therefore a 
polymer chain is expected to collapse to a small region 
in space. However the presence of excluded volume of 
the chain prevents this from occurring. The statistics of 
a chain in this case was the subject much controversy. 
This was settled by the work of Cates and Ball [l[ who 
showed that the statistics of a chain on an infinite lattice 
with frozen disorder, that is the quenched case, was iden- 
tical to that of a chain with mobile defects, that is the 
annealed case. They showed this by dividing up the infi- 
nite lattice into large finite regions with frozen disorder. 
The average properties of a chain are obtained by averag- 
ing over all these large finite regions. This is equivalent 
to performing an annealed average. It is easy to perform 
the annealed average over uncorrelated random disorder 
and show that this has no effect of chain statistics. Sim- 
ulations Q away from the percolation threshold bear out 
this prediction. Work on off-lattice models of spherical 
obstacles show more complicated behavior, because the 
correlations in the disorder are no longer point-like, so 
this influences the conformation of a chain. For low ob- 
stacle volume fraction, these correlations induce an ef- 
fective attraction causing a decrease in chain size. As 
the volume fraction gets close to the percolation thresh- 
old [1, Q , the exclusion of phase space for finite lattice 
sizes leads to an increase in the average chain size. 

Although chain statistics in the above situation are 
not affected by frozen disorder, there is another impor- 
tant physical quantity that does change. Consider the 



density of monomers for a given background of disor- 
der. If this is averaged over time, one expects that the 
density will fluctuate because of the fluctuations in the 
random potential. This problem was studied numeri- 
cally by Gersappe et al Q in two dimensions Q where it 
was shown this time averaged monomer density (TAMD) 
obeyed multifractal statistics fl2l - fl4j over the parameter 
range studied. The spectrum of multifractal dimensions 
was computed numerically and gave strong evidence for 
multifractal scaling. 

More recently [15] , the properties of a polymer chain on 
the backbone of a percolating cluster were examined nu- 
merically at the percolation transition. They found that 
it exhibited multifractal properties in agreement with the 
theoretical work of Janssen and Stenull @ . A measure of 
bond density for an ensemble of chains of varying lengths 
that connect two points separated by a distance R was 
calculated numerically. That measure was analyzed and 
shown to have a spectrum of multifractal exponents. 

The earlier work of Gersappe et al Q in two dimen- 
sions claimed that multifractal statistics were true even 
far from the percolation transition where the disorder is 
dilute. In other words, the multifractal nature of poly- 
mers in disordered systems is much more general and in 
fact, a simpler more experimentally accessible quantity, 
namely the time averaged monomer density, can be used 
as a measure. 

Multifractal distributions are characterized by very 
large fluctuations, so in this case, we expect that the 
TAMD p(r) will have a probability distribution that be- 
comes increasingly broad as the system size in increased. 
We also expect multifractal statistics to be present for 
scales less than the average radius of gyration, R g , of 
the chain. On larger scales, the density fluctuations 
should saturate so that the distribution P of p should 
obey (HQ 

P(logp) oc exp(log(i? 5 )/(log(p)/log(i? s )). (1) 

In this paper, we study this problem in two dimen- 
sions for different amounts of disorder and show that 
multifractality is seen in all cases. The original work 
of Gersappe et al only presented results for one value of 
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disorder. Next we study this system in three dimensions 
and show that multifractal statistics persist for this case 
as well. Next we consider many chains at finite concen- 
trations and demonstrate how this smooths out density 
fluctuations so that the statistics are no longer multifrac- 
tal. 

In section |TT] we describe the model that we will apply 
to study this problem. In section HIIl we study the TAMD 
in two dimensions. The results are similar to what was 
found earlier [3] but with much improved data providing 
more solid evidence for mult ifr act ality of a range of ob- 
stacle densities. Section HVl considers the same problem 
in three dimensions. Although box sizes are smaller, this 
work provides solid support for the case the TAMD is 
also a multifractal measure in this case as well. In sec- 
tion [V] this problem is analyzed for a solution of chains 
at finite concentration where it is clear that the system 
ceases to be multifractal over the scale of the radius of 
gyration of a chain. An interesting bimodal distribution 
for the distribution of the TAMD is found. 



that length scale, then it is straightforward to show that 

(L/l) d 

T q {i) = C£ #>~rw-* (3) 

r 

By analyzing how this sum depends on I, r{q) can then 
be determined. 



III. SINGLE CHAINS IN TWO DIMENSIONS 

The above model was studied for a single chain in two 
dimensions for different values of the disorder and chain 
lengths. The TAMD p r is plotted in Figure Q] Figure 
[IJa) shows a chain of length TV = 512 in a box of length 
L = 128 with an obstacle density of n — .1. The self 
similar nature of the measure is apparent. Figure [Tfb) 
shows the TAMD for N = 256 and L = 64, with n a = 0.2. 



II. THE MODEL 

We consider a cubic lattice in d = 2 and 3 dimen- 
sions, with m self avoiding walks of unit step length of N 
monomers. The monomers move in a cubic box of width 
L. Skew boundary conditions are employed [l6j]. Ob- 
stacles are randomly placed on sites which are excluded 
from occupation by polymer chains. We denote the av- 
erage fraction of excluded sites as n Q . 

The chains are moved using "reptation dynamics" [l7j | . 
Briefly, the head or chain of a chain is picked at random 
and an attempt is made to move it to the other end. 
The attempt fails if the site is occupied by a chain or 
an obstacle, otherwise it succeeds. After equilibrating, 
statistical properties of the chains are measured, such as 
the radius of gyration and the time averaged monomer 
density. 

To test this method, we checked that the internal chain 
statistics of a single chain with random disorder were 
identical within statistical error to that of a self avoiding 
walk with no disorder. 

Next we considered density fluctuations. Moments of 
the normalized time averaged monomer density p r can 
be calculated by summing moments of the density over 
all lattice sites and inferring how this depends on lattice 
size 

L d 

(£p q r)~L-^ (2) 

r 

The average is over different realizations of the disorder. 
However, it is more accurate to obtain the exponents r(q) 
by coarse graining, making use of the self-similar nature 
of multifractals. If we coarse grain over boxes of width I, 
and define p r as the density measure coarse grained over 




(b) 

FIG. 1: (Color Online) The time averaged monomer density 
for a self avoiding random walk in two dimensions, (a) N = 
512, L = 128, and n = 0.1; (b) TV = 256, L = 64, and 
n = 0.2. 



The density shows a self similar structure with different 
regions appearing statistically quite similar but differing 
by an overall multiplicative constant as expected for a 
multifractal measure. For distances larger than the ra- 
dius of gyration, one expects multifractal scaling to no 
longer hold. Similarly, for distances smaller than some 
cut-off distance Iq, related to the density of obstacles, the 
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density should become smooth so that multifractal scal- 
ing no longer applies. But over intermediate scales, as 
can be seen, the density appears self similar. This will 
be now quantified below. 

Coarse graining of p is done by recursive decimation of 
the original lattice into 2x2 blocks. This implies that the 
coarse graining length is a power of 2, I — 2 m where m is 
the number of times the lattice is decimated. Good fits to 
power laws were obtained between the lower and upper 
cutoffs for Eq. [3] as shown in Fig. [2] This plots the left 
hand side of Eq. [3] as a function of coarse graining m = 
log 2 l- Fig. [2]shows the case N = 512, L = 128, and n = 
0. 2 for q = — 1,0,1,2. The results were averaged over 200 
different realizations of the random obstacles, and each 
simulation was run for a total of 8 x 10 9 steps. This is 
strong evidence of self similar multifractal behavior far 
away from the percolation threshold in two dimensions. 
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FIG. 2: (Color Online) The scaling of moments of the coarse 
grained density T q as defined in Eq. [3] as a function of the 
amount of coarse graining m in two dimensions. The top line 
is data for q = — 1, the middle, q = 2, and the bottom is for 
9=1- 



Using these hts, r(q) defined in Eq. [3] was calculated. 
The results are shown in Fig. [3] For clarity we plot the 
difference in r defined as At (q) = r(q) — r(q— 1/2). This 
is shown for three different sets of parameters, no = 0.1, 
0.2, and 0.3. 

The data give strong evidence for multifractal behavior 
for a range of different obstacle densities. 



IV. SINGLE CHAINS IN THREE DIMENSIONS 

The simulation was run in three dimensions and ana- 
lyzed for the possibility of multifractal scaling, following 
the same procedure as in two dimensions. We first plot 
(p q ) as a function of the coarse grained length scale I, as 
is defined in relation to Eq. [3j Fig. [4] shows the case 
N = 128, L = 32, and n — 0.3. Because the box size 
and the chain length are smaller than in two dimensions, 
power law scaling can only be seen over a more limited 
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FIG. 3: (Color Online) The results of simulations to de- 
termine the multifractal dimensions for for different obstacle 
densities no = 0.1 (triangles), 0.2 (squares), and 0.3 (dia- 
monds) in two dimensions. 



range, however, it does support the hypothesis that the 
TAMD is multifractal. For q — 2 a slight curvature can 
be seen in the log-log plot which is to be expected from 
finite size effects. Overall the fits to power laws over this 
range is very good. Similarly good fits are also seen for 
no = 0.2. In Fig. 03 we display fits for At as a function 
of q for no = 0.1, 0.2, and 0.3. 
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FIG. 4: (Color Online) The scaling of moments of the coarse 
grained density T q as defined in Eq. [3j as a function of the 
amount of coarse graining m in three dimensions for no = 0.3. 
The top line is data for q — —1, the middle, q = 2, and the 
bottom is for q — 1. 



Therefore is appears that in three dimensions, our 
simulations and analysis also support the idea that the 
TAMD is multifractal for obstacle densities away from 
the percolation. 
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FIG. 5: (Color Online) The results of simulations to de- 
termine the multifractal dimensions for for different obstacle 
densities no = 0.1 (triangles), 0.2 (squares), and 0.3 (dia- 
monds) in three dimensions. 
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FIG. 6: (Color Online) The distribution of the TAMD p for 
40 chains, an obstacle density of no = 0.2 in three dimensions 
with N = 128, L = 32. 



V. MANY CHAINS 

The internal statistics of a self avoiding walk (SAW) 
in a much larger lattice of quenched random obstacles 
was shown to be equivalent to that of an annealed aver- 
age and hence that of a SAW without disorder For 
the case of a many chain solution, extending this result 
is not completely straightforward. The difficulty in this 
case is that the polymer solution extends throughout the 
entire lattice. However, this problem is easily amenable 
to numerical simulation. Various equilibrium internal av- 
erages for chains were calculated and found to be almost 
entirely unaffected by the disorder. As an example, we 
considered m — 40 chains each of length N = 128 with 
an obstacle density of no = 0.2 in a box of size L = 32. 
Excluding obstacles, this corresponds to a chain filling 
fraction of 20%. R 2 g = 49.190 ± 0.013 in this case. In 
comparison R 2 g = 48.785 ± 0.006 when run with no ob- 
stacles. The difference between these two cases is less 
than 1%. As the obstacle density approaches the per- 
colation point, one would expect the difference between 
these two cases to increase due to the increasing likeli- 
hood of islands that are inaccessible, similar to the single 
chain case 0, H| . 

In contrast to internal chain statistics, such as the aver- 
age radius of gyration, we expect that the TAMD for the 
many chain case to be greatly affected by the presence of 
disorder, as it was for the single chain case. 

For a system of chains at finite density in three dimen- 
sions, the fluctuations in chain density will be quite differ- 
ent than for a single chain. Above the correlation length, 
correlations in density will decrease rapidly. Therefore 
the distribution of the TAMD is no longer expected to 
be multifractal over these scales. This problem was sim- 
ulated using the same approach as above. 

The distribution of the TAMD was averaged over lat- 
tice sites and different realizations of the obstacle density 



in three dimensions. This is displayed in Fig. [S] for the 
same system as discussed above: m = 40 chains each of 
length N = 128 with an obstacle density of no = 0.2 in a 
box of size L — 32. The distribution is plotted with log- 
arithmic axes (base e). The results were averaged over 
155 different realizations of the random obstacles, and 
each simulation was run for a total of 8 x 10 9 steps. The 
TAMD p was normalized so that (p) = 1. 

The results show that the distribution P(log p) has two 
peaks. One peak is close to p = 1 and most values of the 
density are clustered around that value. However, there 
is another peak for very small values of the obstacle den- 
sity p = 0.003. The second peak is unexpected because it 
implies that there are intermediate values of the density 
that are much less probable. 

The same feature persists when the obstacle density 
no is reduced to 0.1 and in both cases it is clear that this 
second peak is not due to random statistical fluctuations. 
This is likely due to the presence of some particular local 
configurations of obstacles, for example a cul-de-sac. In 
this case there will be entropic exclusion of chain for such 
a region leading to a lower monomer density. 

The distribution in Fig. [S]is is not compatible with the 
prediction of multifractal statistics (Eq. [T|) as the distri- 
bution must be convex. In addition, the coarse graining 
analysis of the last two sections does not yield a spec- 
trum of exponents. Ar(g) over the same range as Fig. [5] 
varies by less than 0.05. The lack of large fluctuations is 
to be expected due to the presence of a screening length. 
Although it is likely that the TAMD is multifractal below 
this scale, it is not possible to investigate this numerically 
because the accessible sizes only cover a small range in 
scales. 
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